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Ultracold atoms in optical lattices offer a great promise to generate entangled states for scalable 
quantum information processing owing to the inherited long coherence time and controllability over 
a large number of particles. We report on the generation, manipulation and detection of atomic 
spin entanglement in an optical superlattice. Employing a spin-dependent superlattice, atomic spins 
in the left or right sites can be individually addressed and coherently manipulated by microwave 
pulses with near unitary fidelities. Spin entanglement of the two atoms in the double wells of the 
super lattice is generated via dynamical evolution governed by spin superexchange. By observing 
collisional atom loss with in-situ absorption imaging we measure spin correlations of atoms inside 
the double wells and obtain the lower boundary of entanglement fidelity as 0.79=b0.06, and the 
violation of a Bell’s inequality with £= 2.21 ±0.08. The above results represent an essential step 
towards scalable quantum computation with ultracold atoms in optical lattices. 


During the last decades, quantum entanglement, the 
key resource for quantum information processing[I], has 
been created in many systems like photons [2], ions|3], 
superconducting circuits [4], and solid-state qubits [5]. 
They have been wildly used for studying quantum 
computation[g] and quantum simulation 7]. Nowadays, 
a significant demand towards scalable quantum infor¬ 
mation processing is to efficiently construct multipar¬ 
tite entangled states. Since ultracold atoms in an op¬ 
tical lattice [8] have excellent coherence properties and 
can be manipulated in parallel, an attractive protocol0 
was proposed to create resilient entangled states for 
measurement-based quantum computation 10, 11 with 
optical superlattices. In this protocol, maximally entan¬ 
gled Bell-type states are first prepared in double well ar¬ 
rays (DWs) of a superlattice, then these Bell pairs are 
connected to each other to create cluster states by using 
Ising-type superexchange interactions|l2], and finally a 
computational algorithm is implemented by performing 
single-particle measurements together with unitary oper¬ 
ations. 

Along this direction, a preliminary step has been taken 
in an exciting experiment by observing and controlling 
the superexchange interactions between two atomic spins 
in optical superlatticesp2]. However, due to difficulties 
for addressing single spins and measuring spin correla¬ 
tions in superlattices, it remains a challenge to generate, 
characterize, and manipulate the entanglement of atomic 
spins, which are the essential ingredients for measure¬ 
ment based quantum computation [8j |9]. 

In this work, we developed a new superlattice con¬ 
figuration featured with spin dependence, which allows 
us to individually address and manipulate the spins in 


the left or right sites of the DWs with microwave (MW) 
pulses. Relying on this configuration, high-fidelity ini¬ 
tialization of the DWs from ||,|) to |t, i) is achieved 
(|t) and | \) denote the two spin states while the comma 
separates the left and right occupations). The entangled 
state, (|t, 4-) + *|4-,T))/v / 2, is generated via a VSWAP 
operational [14] . Afterwards, an effective short-range 
gradient magnetic field is induced to manipulate the 
phase of the entangled state and transfer it to the Bell 
state (| t> 1} ± ||, f))/\/2. Then, employing a new de¬ 
tection routine with two-stage filtering and imaging, we 
measure quantitatively the spin correlations, (S z <g>S z ) 
and (Sy®S y ) of the two atoms by observing the atom loss 
arising from hyperfine changing collisions [15j 16]. The 
lower boundary of the entanglement fidelity and the vi¬ 
olation of the Clauser-Horne-Shimony-Holt (CHSH) Bell 
inequality [171 are finally derived from the spin correla¬ 
tions. 

Here we stress that the atom cloud is prepared in a 
pancake-shaped trap deeply in 2D regime on the X-Y 
plane, therefore the longitudinal inhomogeneity along the 
Z direction existing in 3D lattices is excluded. Moreover, 
an in-situ absorption imaging along the Z direction is em¬ 
ployed for detecting the density distribution of the atoms, 
which has a great advantage of suppressing the transver¬ 
sal inhomogeneity (in X-Y plane) by selecting a proper 
area in the sample. The two approaches are essential for 
generating entangled states with high fidelity and long 
coherence time. The present techniques and experimen¬ 
tal routines form an architecture for generating scalable 
cluster states [9] as well as many-body Hamiltonians [T8], 
therefore are ideally suited for studying measurement- 
based quantum computation [10] and performing quan- 
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FIG. 1. The experimental apparatus, (a) Three retro-reflected optical lattices create the double well arrays: a short lattice 
V XSJ a long lattice V x i along the X direction, and lattice V y along the Y direction. The imaging beam is applied along the Z 
direction, (b) The laser beams around the atom cloud. The 2D gas is loaded into a single layer of a pancake-shaped trap, 
which is a series of 4/im-distance layers created by interfering two beams, (c) The density distribution of the MI derived from 
an average over 10 samples by in-situ absorption imaging. The central part of the atom cloud with a diameter of 16 pixels 
(0.93 /im per pixel) is used as the region of interest (ROI) for later studying spin dynamics and entanglement. The averaged 
filling in the ROI is about 0.8 atoms per site. The geometries of the lattices are shown on the X and Y direction. 


turn simulation [T9]. 

The system under consideration is a series of isolated 
double-well arrays created by an optical superlattice — 
a system which can be well described with the Bose- 
Hubbard model characterized by a nearest-neighbor tun¬ 
neling J, onsite interaction energy U and the imbalance 
offset A in DWs[T2], All these parameters can be well 
controlled by the intensities and frequencies of the lattice 
lasers. In the following part we focus on the condition of 
balanced DWs (A = 0) and deep lattices where the inter¬ 
action dominates U J. With the initial states in the 
subspace of singly filled states of | T, 1} and | V T), double 
occupation at each site is allowed only virtually. Then 
the Bose-Hubbard model can be described as a Heisen¬ 
berg spin model H = — J ex ^L -S#, where J ex =4:J 2 /U de¬ 
notes the superexchange coupling between the two sites 
in a DW. By initializing the DW to | T, 1} and using ap¬ 


propriate potentials, one can observe the superexchange 
driven evolution from |T,i) to |^,T), and vice versa[T3]. 
A VS WAP operation is applied to generate the entangled 
state V) = (|tV)+^|V t))/y/2 by halting the superex¬ 
change driven evolution at ti = h/AJ exi with h being the 
plank constant. However, we intend to prepare the Bell 
state |t) = (|t4) + |l,t))/v / 2 (the spin triplet), which is 
taken as the starting point for creating cluster states for 
measurement based quantum computation|j]. In order 
to achieve |£) from V), an additional phase between the 
two components |T, \) and |i, t) has to be introduced. 
We first induce an effective gradient magnetic field in¬ 
side the double-well by creating a spin-dependent [14] su¬ 
perlattice. This causes non-degeneracy between the two 
components with an energy split of 5. By holding the 
atoms in this lattice for a period £ 2 , an additional phase 
of S -£2 is accumulated. Therefore the phase of the en- 
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FIG. 2. Observation of spin dynamics driven by the superexchange interaction, (a) Initialization, evolution and 
detection of spin states in the spin-dependent superlattice. The lattice potentials for |^) and |t) are different, and the coupling 
frequency of |>|r) *<=> |t) i n the left site is shifted 31.8 kHz away from that in the right site. The spins in left sites can be 
individually addressed with a MW pulse, and the DW can be initialized from |^,^} to Then the spin state is evolved to 

a superposition state by superexchange driven evolution. By flipping the spin in every left site before the imaging pulse, the 
DW states of \7,l) and ||,t) are transferred to and respectively, for detection. Only the |t) can be coupled with 

the imaging pulse. With this detection method, the superexchange driven evolution is observed by recording the population 
of 11) in two experimental conditions: (b) V xs =12 Er, V x i =10 Er, J/U = 0.37, a numerical simulation of the BHM agrees 
well with the experiment data; (c) V xs =20 Er, V x i =10 Er, J/U = 0.04, a damped sine fit shows a superexchange frequency 
of 10.0(5) Hz with a 1/e lifetime of 280(20) ms. 


tangled state can be modulated and the spin triplet \t) 
is obtained by controlling t 2 . 

Our experiment starts from preparing a two dimen¬ 
sional (2D) quantum gas by loading a nearly pure 87 Rb 
Bose-Einstein condensate (BEC) into a single layer of a 
pancake-shaped trap, which is created by interfering two 
laser beams with wavelength of A s =767 nm and an an¬ 
gle of intersection of 11°, as shown in Fig.l(a,b). The 
confinement of this trap is highly anisotropic, with trap 
frequencies oo x ,y,z « 27 rx (16,14,7000) Hz, leading to an 
aspect ratio of 467:1. This 2D sample contains around 
1.2 xlO 4 atoms spin-polarized in | \) and reaches a tem¬ 
perature as low as T 2 d = 23(2) nK. We use the hyperfine 
levels 1 1) = \F = 1, mp = —1) and 11) = \F = 2, mp = 
—2) to represent the pseudo spin-1/2 system. The ratio 
hbJz/k^Ttt ~ 14.6 with k b being the Boltzmann constant, 
indicates that the sample is deeply in the 2D regime [20]. 

Then the atoms are adiabatically loaded into a square 
lattice, which is composed of two retro-reflected lattices 
V xs (short lattice) and V y (Y lattice) along the X and Y 
directions with wavelength X xs = X y = 767 nm, as shown 
in Fig.l(a,b). The system enters the Mott insulator 
(MI) regime [2TH24] where most sites are filled with sin¬ 
gle atoms when the lattices are ramped to V xs = V y = 26 
Er, with Er = h 2 /2mX 2 s being the recoil energy and m 
the mass of the atom. Afterwards the two lattices are 
increased further to V xs = 60 Er and V y = 40 Er to freeze 


out the atom tunneling. The density distribution of the 
sample is obtained by in-situ absorption imaging along 
the Z direction with a microscope objective (N.A. = 0.48) 
and a low noise CCD camera, shown in Fig. 1(c). An 
average filling number of 0.80 atoms per site in the cen¬ 
ter part of the atom cloud is obtained. Measuring the 
atom loss of hyperfine changing collisions in the lattices 
reveals that 76% of the lattice sites are filled with one 
atom, about 2% are filled with two atoms and the rest 
are vacant sites [25]. A great advantage of using in-situ 
imaging is that one can select a proper area in the sam¬ 
ple with optimal filling properties and least transversal 
inhomogeneity (see later text). 

Next, the MI atoms are transferred into the DWs of 
a superlattice along the X direction. By superimpos¬ 
ing another retro-reflected lattice V x 1 (long lattice, wave¬ 
length X x i = 1534 nm) with the short lattice, a su¬ 
perlattice V x (x) = V XS cos 2 (k x x) — V x i cos 2 (k x x/2 + ip) is 
formed, where k x = 2tt/X xs is the wave vector and ip is the 
superlattice phase. Controlling the relative frequency be¬ 
tween the long and short lattices, the superlattice phase 
is tuned to 0, therefore balanced DWs with spin state 
| j,, j.) are prepared. 

We then initialize all the DWs to 11, i) by flipping the 
atomic spins in every left site inside the spin-dependent 
superlattice with a MW pulse. Concretely, tuning the 
voltage applied to the electro-optical modulator (EOM), 
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FIG. 3. Transversal inhomogeneity observed in superexchange evolution, (a) Experimental time evolution of line 
density profiles during the superexchange dynamics in the condition of V x i = 10 Er and V xs = 16 Er (J/U = 0.11). (b) Each line 
density profile is the mean cross section along the Y direction by averaging the rectangular region (14.9 /xm width) along the 
X direction. A dephasing process is clearly visible along the Y direction, (c) The simulation of experiment with the measured 
laser diameters w xs = 240 /xm, w x i = 300 /xm and w y — 240 /xm. The evolution of the superexchange resolved by using either the 
whole atom cloud (d) or using the center part (diameter of 14.9 /xm) (e) as the ROI for counting atoms. In the first case, a 
damped sine fit gives a 1/e life time of 41(5) ms, while 115(12) ms for the second one. 


an angle of 0 between the polarizations of the incident 
and the retro-reflected beams of the short lattice is cre¬ 
ated, therefore a spin-dependent superlattice is built 
up [25]. In the condition of V x i =56.3 Er, V xs = 150 Er and 
6 = 46°, the resonant microwave frequency for coupling 
I 1}l O | t)z, in the left site is shifted 31.8 kHz away from 
the resonance | \)r | | )r of the right site, shown in 

Fig.2(a). A high addressing fidelity is achieved by using 
a MW 7r-pulse with Rabi frequency of D = 2tt x 8.1 kHz, 
during which the magnetic field B x = 95 /xT along the X 
direction is actively stabilized and the noise is suppressed 
to less than 16 nT. 

After state initialization, the spin dependence is 
switched off by ramping down the EOM voltage. Mean¬ 
while the quantum axis is rotated from X to Z direction 
to avoid the residual magnetic gradient along the X di¬ 
rection. The superexchange is now initialized by first 
ramping down the long lattice to 10 Er, and afterwards 
ramping down the short lattice to the final value Vse in 
500 fjbs. After letting the system evolve for a time G, we 
halt the superexchange by ramping up the short lattice 
to 60 Er in 500 /xs. The spin configuration is then freezed 
until the state detection. 

To detect the spin state inside the DWs, the quan¬ 
tum axis is aligned to the X direction, then the spin- 
dependence is switched on for addressing the left sites. 
Flipping the spins in the left sites, the two spin con¬ 
figurations will end up along two different channels: 
|t,l) => ||,|) and ||, t) => as shown in Fig.2(a). 

Only the 11) state can absorb the imaging light (cycling 
transition \5Si/ 2 ,F = 2) <t=> |5P 3 / 2 ,F = 3)). Therefore, 
the superexchange-driven oscillation between 1|) and 
11, |) can then be observed with the in-situ imaging by 


comparing the population of atoms in spin up to the 
total atom number in the ROI N t . Two typical oscilla¬ 
tions measured at different DW potentials are shown in 
Fig.2. For a low barrier height (J/U = 0.37), the dy¬ 
namics can be directly described with the Bose-Hubbard 
model, in which the fast oscillations correspond to the 
first-order tunnelings. The simulation matches well with 
the experiment as shown in Fig.2(b). For larger barrier 
heights (J/U <C 1), the oscillation is then mainly driven 
by the superexchange interaction. In the experimental 
condition of V x i = 10 Er, V X3 = 20 Er (J/U = 0.04), an 
oscillation frequency 10.0(5) Hz and a decay constant of 
280(20) ms is obtained. 

Here we study the dephasing effect of the superex¬ 
change evolution. As shown in Fig.3(a,b), the dephas¬ 
ing along the Y direction is clearly observed when we 
perform the superexchange in the DWs with V x i = 10 
Er and V xs = 16 Er. A simulation of the Bose-Hubbard 
model with considering the lattice profiles is shown in 
Fig.3(c), which agrees well with the experimental data. 
We then find that the dephasing process is mainly caused 
by the transversal inhomogeneity of the lattice profiles 
through spatially modulating the Bose-Hubbard param¬ 
eters of the DWs. By using the in-situ imaging system, 
one can resolve the atom distribution with a resolution 
around 2 /xm. Thus we can choose an area in the atom 
cloud as the ROI with least inhomogeneity to obtain a 
long coherence time while keep a good signal to noise ra¬ 
tio. Here, an area with diameter 14.9 /xm in the center of 
the atom cloud is used as the ROI. The resulting coher¬ 
ence time of the superexchange is 3 times longer when 
using the ROI instead of the whole atom cloud, as shown 
Fig.3(d,e). The heating effect due to scattering the 767- 
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nm lights of the pancake and lattices is considered as 
another limitation of the coherence. 

The VS WAP operation is realized by halting the 
superexchange-driven evolution (V x i = 10 Er, V xs = 20 
Er) at t\ = 25.1 ms, and the DWs are then prepared in 
the entangled state | ip). In order to control the phase 
between the two components in this entangled state, we 
induce an effective short-range gradient magnetic held 
inside the DW by switching on the spin dependence with 
V x i = 5.6 Er, 14 s = 60 Er and 0 = 7.5°. This gradient cre¬ 
ates a non-degeneracy of £ = 27rx427 Hz between the 
two components and initiates an oscillation 26] between 
the spin singlet state (|s) = (|t? 1) ~ \iA ))/and the 
triplet state 1£). After a singlet-triplet oscillation (STO) 
time £ 2 , we switch off the spin dependence and ramp 
down the long lattice to 0. By measuring in the rotated 
basis |±) = (|t) ±i|4'))/v / 2, i-e. applying a 7r/2-pulse on 
both DW sites, the triplet and the singlet state can be dis¬ 
criminated: 1 1) (|t, t) + | V t))/V2 and \s) | s). 

Then a two-stage hltering and imaging routine is em¬ 
ployed to determine the entanglement phase, as shown in 
Fig.4(a). We assume the complete spin correlated basis 
in the DW, {||,|), |t4), IVt), First, an imag¬ 

ing pulse is used to count ah the spin-up atoms and 
remove them from the lattices by heating. Afterwards, 
a MW rapid adiabatic passage is applied to hip ah the 
remaining atoms to |V) = \F = 2, mp = —1)- Then, the 
atoms are transferred to the long lattice of 35 Er in 2 ms 
by ramping down the short lattice. Meanwhile the lat¬ 
tice depth along the Y direction is increased to 60 Er for 
enhancing the onsite interactions. With this condition, 
the double occupancies will be removed from the lattice 
after a holding time of 500 ms due to hyperhne chang¬ 
ing collisions, while the single fillings will survive. Then 
the in-situ imaging is performed for counting the remain¬ 
ing number of atoms AY By this method, the triplet 
state can not contribute any count to the final number of 
atoms, while one atom of the singlet will be counted. As 
shown in Fig.4(b), an oscillation curve with a period of 
2.34(2) ms is obtained by counting the number of atoms 
in the ROI at different STO times with this detection 
routine. According to the STO curve, the phase of the 
entangled state can be well controlled and the Bell state 
1£) can be prepared after a STO time £2 = 0.3 ms. 

Following the two-stage hltering and imaging routine, 
we verify the spin-entanglement in the DWs by com¬ 
paring the sum of two images (Ah + £V 2 ) to the to¬ 
tal number of atoms N t in the ROI. The number of 
atoms in the correlated basis of | V is then derived 
as = N t — (Ah + N 2 ). By addressing and hipping 
the spins on the left site, right site or both sites before 
the detection process, we can project the states | j"V), 

| Vt)> or | t?t) the measurement state |t,t) and de¬ 
rive the number of atoms for £V^, N^, or Ahq, re¬ 
spectively. From these measured spin fractions, we can 
get coincidence-like probabilities P a ^ = N a ^/ v N u , v 

(a,b,u,v or and the spin correlation of ( S z ® S z ) 
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FIG. 4. Measurement of the entanglement phase, (a) 

Two-stage hltering and imaging routine for detecting the 
phase of entanglement: i) the complete spin correlated ba¬ 
sis in the DW; ii) remaining | V) atoms after recording the 
number of | f) atoms and the MW rapid adiabatic passage; 
iii) merging the DWs and holding the atoms for 500ms in a 
deep lattice to remove the double occupancies; iv) counting 
the number of remaining atoms AY The four different ini¬ 
tial pair states will contribute differently to the hnal atom 
counting, thus by projecting the entangled state to |±) basis 
one can resolve the singlet and triplet state, (b) Time evo¬ 
lution of the remaining atom number for different STO time 
after the two-stage hltering and imaging process. A period of 
2.34(2) ms is derived from a sine ht. 


is E zz = P ^ - P Y - Ppt + Pj^q. 

This technique for measuring spin correlations is sim¬ 
ilar to the measurement of polarization correlations in 
photonic entanglement [2]. Similarly, the spin fractions 
in the |+)/| —) basis can be measured by applying a 
7r/2-pulse to both sites before these measurements, and 
the probabilities P c ^ (c,d = + or —) are derived. The 
spin correlation of (S y ® S y ) is E yy = — P+,+ + P +? _ + 
P_ 5+ — P-,-. From these measurements, as shown in 
Fig.5(a,b), we obtain the lower boundary of the entan¬ 
glement hdelity|27] P > ~{E ZZ + E yy )/ 2 = 0.79 d= 0.06, 
higher than the classical limit of 0.5 by 5 standard devi¬ 
ations. 

Furthermore, we demonstrate the existence of spin en¬ 
tanglement by violating the CHSH-type Bell inequality. 
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FIG. 5. Measurement of spin correlations and entanglement, (a-b) The coincidence-like probabilities derived in 
|t>/|;> and |+)/|-> basis, P t>t = 0.05(3), P t ,t = 0.45(3), P ut = 0.45(3), and P Ui = 0.04(3); P+,+ = 0.44(5), P+,_ = 0.06(4), 
P_,+ = 0.05(4), and P-,~ = 0.44(5). (c) Measured spin correlation curves of E zz (/3) and E yy (/3). We derive the violation of 
the CHSH type Bell inequality as 2.21 d= 0.08 from the sine fittings. 


The spin-correlation curves are measured by first assign¬ 
ing the left site an additional rotation (3 = <a7r/4 (a = 
1,2, 3..., 8), and then repeating the above measurements 
for the spin correlations of E zz (f3) and E yy (/3), as shown 
in Fig.5(c). Choosing the rotation of <a7r/4 for the left site 
causes minimum response to the atoms in the right site. 
From fitting the two spin-correlation curves, the quan¬ 
tity S = \E(0 u 0 2 ) ± E(0 u 0' 2 ) - E(0[,0 2 ) ± E(010' 2 ) I, 
with (0i, 62 , 0' 2 ) = (0, 7t/2, 37t/4, 57t/4), is derived as 

S = 2.21 ±0.08, which violates the CHSH inequality 
with 2.7 standard deviations. 

The observed entanglement fidelity is limited by sev¬ 
eral experimental imperfections. First, due to the cou¬ 
pling with the environment, the spin states experience de¬ 
coherence during the period of generating the entangled 
states and another holding time of stabilizing the mag¬ 
netic field for further microwave operations. This deco¬ 
herence causes a degradation of the fidelity by about 17% 
. Second, the errors of microwave operations decrease the 
fidelity by about 3%. Therefore the expected fidelity of 
80% is consistent with the experimental measurement. 
Although the generation of spin entanglement in dou¬ 
ble wells is deterministic, extending to longer chains of 
spin entanglement suffers from lattice defects due to finite 
temperature. Further cooling[28] down the temperature 
of the atom cloud will be helpful for suppressing vacant 
sites and achieving unit filling in the lattices. 

In summary, we have demonstrated the generation, 
manipulation and detection of atomic spin entanglement 
in an optical superlattice, the first step towards measure¬ 
ment based quantum computation. The spin dependence 


built-in superlattice brings great flexibility for addressing 
the single spins in the double-well arrays, which leads 
to high-fidelity state initialization and detection. Both 
the longitudinal and transversal inhomogeneities are well 
suppressed by confining the ultracold atoms in a 2D plane 
and employing in-situ imaging. Therefore, a long coher¬ 
ence time of the entangled state is achieved. The rou¬ 
tine of two-stage filtering and imaging makes it possible 
to measure the spin correlations from which we can de¬ 
rive the lower boundary of the entanglement fidelity as 
0.79±0.06 and the violation of the Bell’s inequality with 
S = 2.21 ±0.08. This method may be used for charac¬ 
terizing large scale entangled states by investigating the 
entropic inequalities[29]. 

By connecting the Bell pairs, one can extend the en¬ 
tanglement to a lattice chain, and further to a 2D plane 
with the presently developed techniques and experimen¬ 
tal routines. Large 2D cluster states can be generated as 
universal resources for quantum computation [9f 3J0]. To¬ 
gether with high-resolution single-site [31, 32 and single¬ 
spin manipulation [33], it may lead to one-way quantum 
computation[lOl |TT] . Owing to the site-resolved spin ad¬ 
dressability in the spin-dependent superlattice, abundant 
physics of many-body systems can be studied with the 
present setup, e.g. ring exchange in a spin plaauette[34] 
of the minimum instance for the Kitaev model [35] and 
quantum magnetism of spin systems [36]. 

During preparation of the manuscript, we became 
aware of a recent related work by T. Fukuhara et al. m\ 
detecting a spin-entanglement wave in a Bose-Hubbard 
chain. 
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SUPPLEMENTAL MATERIALS 


Preparing balanced DWs 

The superlattice potential along the X direction is V x (x) = V XS cos 2 (k x x) — V x i cos 2 (k x x/2 + </?), when the EOM 
modulation is off. By beating the up-converted light of the 1534-nm laser with the 767-nm laser, the relative frequency 
of the two lasers are locked at about 5.5 GHz with a linewidth less than 400 kHz in the frequency offset locking scheme. 
Therefore the relative phase c p between the long lattice and the short lattice can be continuously adjusted from 0 to 
7r by tuning the offset locking frequency. During the measurement of the superexchange, STO and entanglement, the 
superlattice phase is tuned to 0. 


Spin-dependent superlattices and site-resolved spin manipulation 


The spin-dependent superlattice is realized by overlapping a spin-dependent short lattice with a normal long lattice, 
as shown in the setup, Fig. 1(a) in the matin text. The spin dependent short lattice is created by modulating the 
polarizations of the forward and the retro-reflected beam. The initial polarization of the forward beam is linear 
polarized in the horizontal plane (X-Z plane), and it is modulated by the control voltage applied to the electro-optical 
modulator (EOM). The axis of the EOM is aligned to the horizontal plane with an angle of 45° and a quarter-wave 
plate (QWP) with its fast axis in the horizontal plane is used to adjust the circular polarization induced by the EOM. 
Then the beam after the QWP is still linear polarized, but rotated by an angle of 6/2 with respect to the original 
polarization, where 6 is the modulated phase by the EOM voltage. The two forward beams of short lattice and long 
lattice (horizontal polarized) are overlapped on a dichroic mirror and focused to the atoms by an achromatic lens. 
For the retro-reflected beams, a special designed QWP is used, which acts as a QWP for the short lattice while 
inducing no phase change for the long lattice. The fast axis of this QWP is also aligned in the horizontal plane. 
The retro-reflected beam of the short lattice is linear polarized with an angle of — 6/2 to the horizontal plane, and 
—0 to the forward beam (Fig. SI (a)), while the retro beam of the long lattice is still horizontal polarized. Thus the 
interference of the short lattice can be modulated, and the atoms in the spin down (\\) & 5Si/2 \F = 1 ,mf = — 1)) 
and spin up (|t) = 5 Si/2 \F = 2,mf = —2)) state will respond differently due to the different dipole couplings to the 
left/right circular polarization components of the lattice. The spin-dependent term is proportional to i(E*xE)-/x, 
where E is the direction of electric filed and /i is the magnetic momentum [SI]. Therefore the largest spin-dependent 
effect is present when the quantum axis is aligned along the X direction. The modulated lattices along X direction 
for the two spin states can be described as 


V a (x,0)=V X8lO 


■At cos 2 (k x x + + A a cos 2 (k x x - 


~V xit(T cos 2 (^x + ip), 


(SI) 


where a = |4) or \\) and A v a is the response for spin a to the circular polarization component of v (+/— for right/left 
circular polarization). Here we have A+ = 0.55, A^ = 0.45, A^ = 0.40, and A jT = 0.60, respectively. The effective 
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FIG. SI. Spin dependent short lattice, (a) Schematics of the linear polarized forward and retro reflected beams with 
polarization difference of 0 , which is modulated by an EOM and QWPs; (b) Effective short lattice depth varied by the EOM 
modulation for a 60 Er short lattice. The calibration is done by applying a short time parametric modulation to the lattice in¬ 
tensity when the long lattice is off. This confirms the theoretical predications (solid lines) of V^(0) — V XS}(T \/ cos 2 6 + B a sin 2 0, 
(7=|4) (blue) or \l) (red). 
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FIG. S2. Site-resolved MW addressing spectrum in the spin-dependent superlattice. Atoms are loaded in lattices 
V x i = 56.3 Er, V xs — 150 Er and V y = 40 Er for transition of |t) => |>f)- (a) When the EOM modulation is off, the spectrum 
has a normal single maximum peak, which means the left and right sites have the same transition frequency, (b) When the 
EOM is modulated to an angle of 46°, two distinct peaks with 31.8 kHz separation are observed, corresponding to the two 
microwave frequencies for addressing the left or the right sites. (c,d) The DW potentials are shown when the EOM modulation 
induces an angle of 0 — 0° or 0 m 46° between the polarization of the forward and retro-reflected beam in the short lattice. 


short lattice depth will be modulated by the EOM, and it is calibrated by the resonance frequency of the coupling 
from the ground band to the second excited band. The effective lattice depth for | \) and 11) are depending on the 
EOM modulating: (#) = V X8j(T V cos 2 9 + B a sin 2 6 , as shown in Fig.Sl(b). In our experiment, B± = 1.13 x 10 -2 

and B^ = 4.54 x 10 -2 . 

We find that the superlattice phase will also be modulated by the EOM, since V a (x,0) can be written as 
V a (x, 6) =V XS ^ cos 6 cos 2 (k x x) + [l - cos# + (A+ - A~) sin#] 

K » (S2) 

- Vxi,a cos 2 (-^-x + ip) - V xs ,a (A+ - A a ) sin 9 cos 2 (k x x - -) 


The last term in Eq.2 is the part which will introduce the spin-dependence, and it will also change the superlattice 
phase ip. For simplicity, we consider the balanced DWs ip = 0 and an EOM modulation of 0 < 0 < 90°. Since 
A + > , the DWs of 1 1) will be tilted to the right side, while < A^ , the DWs of 11) will be tilted to the left 

side, as shown in Fig.S2(c,d). In this case, the energy shifts for 11)/| 1} in left/right sites of the DWs are different. 
With this method, we can address the left or right sites in the DWs without affecting the other ones. 

Shown in Fig.S2, two distinct micro-wave spectrums are obtain by a MW pulse with Rabi frequency = 3 kHz 
and pulse length 167 /is in the case with or without the modulation of EOM. The spectroscopies are measured at 
lattice depths of V x i = 56.3 Er, V xs = 150 Er by recording the number of atoms flipped by the MW. We first measure 
without EOM modulation, 0 = 0°, and a single resonance peak can be obtained. While two separated peaks with 
separation of 31.8 kHz are observed when we turn the EOM modulation to 9 = 46°. The two peaks correspond to 
the different resonance frequencies for the left and right sites, respectively. 


Measure the filling parameters in MI 

After entereing the MI regime, we increase lattice depth on the X and Y direction to freeze the tunneling. Afterwards, 
we use a rapid adiabatic passage to transfer all the atoms from the initial state ||) to \ Y) (55i/ 2 |F = 2 ,imf = —1)). 
The atoms in the F = 2 mainfold can undergo hyperflne changing collisions (Ref.[15] and [16] in the main text), 
and have a shorter lifetime for the sites filled with two or more atoms. The double fillings will be removed from 
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FIG. S3. Calibration of the filling parameters: the schematics, in-situ distribution and the mean cross section along the 
X direction for the rectangular area (gray shading) with a width of 16 pixels (14.9 /im) for (a,b,c) the initial MI state; (d,e,f) 
the atom cloud after holding 500 ms in the short lattice; and (g,h,i) the atom cloud after the second filtering killing of 500 ms 
in the long lattice. No obvious loss is observed for the first filtering process, which indicates a low probability for filling with 
two atoms, while a clear loss is observed for the second filtering process, which corresponds to the filling of one atom in both 
sites of the DW. 


the trap by this filtering process after a holding time of 500 ms in our experiment. Comparing the atom numbers 
in the ROI with and without collisional loss, we can obtain the double occupancy as j92 ~ 1-8%, as shown in Fig.S3 
(a,b). The average filling in the ROI of the MI is about 0.80, which indicates that multiple occupancies have low 
probabilities. We will neglect the occupancies with more than two atoms for later calibrations and experiments. After 
this filtering, we transfer the remaining atoms to the long lattice by ramping down the short lattice and measure the 
loss after another 500 ms holding time (Fig.S3 (c)). The loss in this process is mainly due to the singly filled states of 
It 7 , tO- Therefore we measured the occupancy of both DW sites filled with one atom as pui ~58%. According to this 
calibration, we know that the filling of 111 in the DWs contain 0.58/0.80 = 72.5% of the total number of atoms in the 
ROI. Furthermore, by assuming that the fillings between the neighboring sites are independent (in deep lattices), we 
derive the filling probability of one atom per site as pi = y/plj7 = 76.2%, and therefore the vacancy filling as po~22%. 


We calibrate the hyperfine changing collisional loss in our experiment by comparing two loss processes: one when 
the lattice sites are mainly filled with one atom and the other one when they are doubly filled. To prepare a sample 
with high single occupancy, we first flip the left sites in the spin-dependent superlattice, remove them from the trap by 
an imaging pulse, and then transfer the remaining atoms into the long lattices (the average filling in the ROI here is 
about 0.4). Then we measure the collisional loss of the atoms in the trap by holding them inside the lattices (V x i = 35 
Er, V y = 60 Er and pancake trap uj z = 7 kHz) for different times. In this case, there is no obvious loss observed since 
there are no collisions for single atom, as shown in Fig.SI3. We intentionally prepare an atom cloud with high double 
occupancy by transferring the MI atoms into the long lattice. The loss process shows an exponential decay with a 
1/e lifetime of 120(5) ms (Fig.S4), and all the double occupancies are removed from the lattice after holding in a 
deep lattice for 500 ms. One may expect a higher collision rate and faster loss process by tightening the confinement, 
however the collision rate in our experiment is mainly limited by the confinement on the Z direction. 
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FIG. S4. Measured collisional loss in the lattices by comparing the atoms in ROI with and without hyperfme spin changing 
collisions: a fast decay with time constant of 120 ms (red) is observed when some of the lattice site are filled with two atoms, 
while almost no decay is observed for the case of only one atom filled in lattice sites (blue), the fitting shows a 1/e lifetime of 
several seconds. 


Calibration of the entanglement phase 


By applying a 7r/2-pulse on both atoms in the DW, the |t)/|l) basis are rotated to the |+)/|—) basis and an 
entangled state \^') = (|t>i) + e l ^\\ n t))/^ with entanglement phase 0 is transferred to the following state as 


W) 




+)) — 


1 - e u 

2V2 




2V2 


e W2 (IT,T> + ll4». 


(S3) 


For the triplet state, (j> is 0, therefore 1 1) (|1bt) + \ While for the singlet state, = 7r, therefore 

IT j 2 

| s) \s). Following the two-stage filtering and imaging routine, the remaining atoms are proportional to the 

singlet fractions, i.e. oc (1 —cos 0). The entanglement phase after a STO modulation is 0(t2) = 0o + ^ 2 , therefore the 
observed atom number for different STO times follows a cosine function in our experiment. 


Calibration of the entanglement fidelity 

The fidelity of the two-atom entangled state generated in our experiment, represented by the density matrix p, with 
respect to the target state 1 1) can be calculated as F=(t\p\t). The lower boundary of the fidelity can be evaluated by 
following the analysis of Ref. [24] in the main text 

F >(Pt,4 +' ^4,4 2 \/^bt ^4,4 + ^+,+ + P-- ~ P+- ~ ^-,+) 

>(^T,4 + ~ -^4,4 + ^+,+ + P-- — — ^-1+) (^4) 

= - (E ZZ + Eyy )/2 


Imaging system and calibration of atomic density 

The imaging objective with an effective numeric aperture of 0.48 contains an aspheric lens and a special designed 
glass plate. The glass plate is used to compensate the wavefront distortion of the glass cell. Another imaging lens with 
a focal length of 500 mm is used to direct the beam into the CCD camera. The amplification factor of the imaging 
system is M = 28.6, and the pixel size on CCD camera corresponds to 0.93x0.93 pm 2 in the real space at the atom 
position. With the absorption imaging of the atoms, the resolution is calibrated to be better than 2.3 pm. 

To faithfully obtain the in situ densities of the atom cloud in the experiment, we apply the saturated absorption 
imaging technique [S2-S4]. A right circular polarized imaging beam with short duration of 10 ps and about 8 saturated 
intensities is applied along the Z direction to the two-dimensional atom cloud. The atomic density of the MI state in 
our experiment is lower than 1 atom per site, which corresponds to an optical density of od~ 2, thus the collective 
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FIG. S5. Calibration of the saturated absorption imaging with the modified Beer-Lambert law by applying different imaging 
intensities to a thermal atom cloud. When a = 2.7, almost the same atom number can be obtained for different imaging 
intensities. 


effect of the atom cloud is small. Furthermore, to avoid multi scattering processes in the highly compressed pancake- 
trap, we release the trap to lower down the 3D atomic density by switching off the lattices and pancake trap for 50 
/is before applying the imaging pulse. 

We calibrate the saturated absorption imaging with the modified Beer-Lambert law S2 


n(x, y)cr 0 = —a In 


h{x,y) 

h(x,y) 


h{x,y) - I t {x,y) 

1 T 

± sat 


(S5) 


where n(x, y ) is the atomic density of the atom cloud, ao is the scattering cross section for the cycling transition of 
5S 1/2 \F = 2) 5P3/2|^ = 3) for circular polarized light, a is a scaling parameter, Ii(x,y) is the intensity distribution 

of the imaging beam, I t (x,y) is the intensity distribution of the imaging beam after being absorbed by the atom 
cloud, and I sat is the saturated absorption intensity. We apply different intensities for the imaging pulse to a thermal 
atom cloud produced in the same condition. By optimizing the parameter <a, the same atom number is obtained for 
different pulses, as shown in Fig.S5. The calibration results a = 2.7(2) for our system. To suppress the large photon 
shot noise at high imaging intensities, a fringe removal algorithm[S5] is applied. 
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